Microscopic theory of photonic band gaps in optical lattices 
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We propose a microscopic model to describe the scattering of light by atoms in optical lattices. 
The model is shown to efficiently capture Bragg scattering, spontaneous emission and photonic band 
gaps. A connection to the transfer matrix formalism is established in the limit of a one-dimensional 
optical lattice, and we find the two theories to yield results in good agreement. The advantage of 
the microscopic model is, however, that it suits better for studies of finite-size and disorder effects. 
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When an electromagnetic wave is sent into an atomic 
cloud, the interference of the radiation fields emitted by 
every atom gives rise to cooperative scattering. In dis- 
ordered systems this interference phenomenon was first 
described by Dicke who evidenced the superradiant emis- 
sion of the cloud [1^ . Later, other striking features linked 
to cooperative scattering have been observed, such as the 
collective Lamb shift [J ^ and a reduction of the radia- 
tion pressure force [5 , 6 . The search for the localization 
of light by the disorder itself is still ongoing [71 [8] . 

The hallmark of the photonic properties of ordered 
atomic ensembles^ such as optical lattices, is the forma- 
tion of a band structure similar to those encountered in 
photonic crystals. Photonic band gaps (PEG) have been 
predicted in one-dimensional arrays of atomic clouds us- 
ing the transfer matrix (TM) technique [9 . This ap- 
proach, which relies on the description of an atomic en- 
semble as a continuous dielectric with a very large trans- 
verse size, describes well the situation of recent experi- 
ments [101412]. which culminated in the first observation 
of a PEG in a one-dimensional optical lattice p!3l [T4]. 

In the case of three-dimensional optical lattices, the 
Eloch-Floquet model has been used to calculate the prop- 
agation of electromagnetic modes in Fourier space and 
identify omnidirectional PEGs in certain geometries as- 
suming infinite and perfectly periodic lattices [13 U^ . 
Nevertheless, omnidirectional PEGs remain to be ob- 
served experimentally. 

In this letter we propose a microscopic model of coop- 
erative scattering from an ordered atomic gas, treating 
the atoms as point-like scatterers interacting with light 
via an internal resonance. We show that this model is 
able to describe the opening of a forbidden photonic band 
due to multiple reflection of light between adjacent lat- 
tice sites. We support our assertion in two ways. Using 
numerical simulations of the microscopic model we find 
that Eragg scattering and PEGs arise in our system. We 
also demonstrate that under a coarse-graining hypothe- 
sis and in the limit of a one-dimensional optical lattice, 
the microscopic model boils down to the TM formalism 
used, e.g., in Refs. [3 HI US]. 

It must be highlighted that our microscopic model does 



not contain the limitations of the above-mentioned other 
techniques. In particular, it does not reduce the atomic 
layers to a smooth dielectric, nor does it assume the 
atomic cloud to be perfectly periodic or infinite in any 
direction. It is thus notably suited to study the role of 
the disorder and finite-size effects on photonic bands. 

The collective light scattering by an atomic ensemble 
is described by the following coupled equations [17[ 118]: 

A A r\ _ ^ . X , r ^ exp(ifeo|r, -rfel) 
V""' - 2 J - 2n^^(^^-) + 2 ^ .^o|r,-r,| 

(1) 

where vj is the position of the jth atom and Pj is the 
excitation amplitude of its dipole. The first term on the 
right-hand side of Eq. ([T]) corresponds to the field Eq (r) 
of the incident laser beam, whereas the last term char- 
acterizes the radiation from all other atoms. Aq is the 
detuning of the incident laser with respect to the atomic 
transition, F is the single atom spontaneous decay rate 
and p is the electric dipole matrix element. 

We test our model on a one-dimensional periodic stack 
of Nd parallel disks randomly filled with Na atoms each, 
illuminated by a laser beam incident under an angle 
with respect to the lattice axis (see Fig.[l]). We compare 
the predictions of the microscopic model with those of the 
TM formalism noting that, while the TM approach as- 
sumes a radially infinite extension of the disks, our model 
is able to account for any distribution, e.g., the Gaussian 
distribution common for thermal atomic clouds. 

We emphasize that Eq. ([T]) describes the scalar light 
scattering. We have also developed a full vectorial model 
giving the results in very good agreement with the scalar 
one derived for the lattice geometry considered in this 
paper. The vectorial model, which will be necessary for 
the description of three-dimensional band structures, will 
be presented elsewhere. 

We first investigate the scattering properties of our sys- 
tem under the Eragg condition, which means that the 
phase-shift of the incident wave between two successive 
atomic disks is tt. In this case, the interference of the 
waves reflected from each disk is constructive, and the 
system is a Eragg reflector. This property is well re- 
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FIG. 1: (color online) Experimental setup: An array of disks 
randomly filled with atoms is irradiated from the side by a 
probe beam incident under an angle Oq. The system can 
be considered as one-dimensional assuming that a,d <C Rd, 
where a and Rd are the thickness and the radius of each disk, 
respectively, whereas d is the inter-layer distance. 



produced by our model in which, despite the point-like 
nature of the scatterers, the incident Gaussian beam is 
reflected by the atomic structure (see Fig. [2|, where we 
consider ^^Rb atoms interacting with the light fields via 
their D2 line. The total electric field E is given by the 
sum of the incident field Eo{r) and the scattered field 
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and according to the extinction theorem, the lattice pro- 
duces in the forward direction a field opposed to the inci- 
dent one. This demonstrates the suitability of our model 
to study, e.g., the microscopic version of the Ewald-Oseen 
theorem [19^ . 

It can be observed that not all incident light is re- 
flected by the atomic structure. A significant part of it is 
re-emitted in the form of the spontaneous emission. This 
phenomenon, which is normally captured in the imag- 
inary part of the refractive index, is naturally present 
in the microscopic model ([T]). The spontaneous emis- 
sion appears in Fig. [3] as the radiation into non-paraxial 
modes. It should be noted that our microscopic model 
does not contain light absorption, and we have verified 
that it conserves energy, i.e., pursuant to the Maxwell's 
equations, the light which is not reflected or transmitted 
is spontaneously scattered more or less isotropically. The 
deviation from perfect isotropy, visible in Fig.|3]as angu- 
lar fluctuations, is a signature of the disorder existing in 
each atomic disk. 

Let us now turn to the study of band gaps. The lat- 
tice reflectivity R= Itat^P and the spontaneous emission 
SE = 1 - R-T (where T = I^atJ^ is the lattice trans- 
missivity) in the microscopic model are in accord with 
the predictions of the TM theory [see Figs. [4]^a), 4(b) 
and Eq.([5|]. Here the reflection and transmission coef- 
ficients are defined as r^^ = Er/Eo and tj^^ = Et/Eo^ 
where Er^t are the total electric fields of the reflected and 
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FIG. 2: (color online) The above picture: Intensity of the light 
in the y — plane as it enters a one-dimensional optical lattice 
and its reflection. The rectangle marks the limit of the atomic 
structure. Below: Zoom of the left part of the atomic lattice. 
The luminous grains correspond to the strong field radiated 
by the atoms close to the y = plane. The simulations are 
realized for N = 9000 atoms randomly distributed over Nd = 
100 layers of thickness a — O.O6A0 and radius R — 9Ao, the 
distance between the atomic disks is d = 0.508Ao with Ao 
being the resonance wavelength. The Gaussian beam of waist 
4.5Ao and power 100 mW is detuned by Ao = F from the 
atomic transition and creates the angle Oq — 0.2 rad with the 
lattice axis. 
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FIG. 3: (color online) Far- field intensity / at the distance 
I5OA0 from the lattice. The light emitted into non-paraxial 
modes exhibits a complex pattern because of the atomic dis- 
order in the disks. The same parameters as for Fig. [2] have 
been used. 



transmitted beam, respectively. Apart from a high reflec- 
tivity, the presence of a PBG is characterized by the local 
density of states (LDOS) vanishing. In the case of a one- 
dimensional optical lattice, the LDOS at the center of the 
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lattice can be conveniently calculated using the complex 
reflection coefficients of the two halves of the lattice 
counting from the lattice center to its ends [20] : 



D = Re 



2 + r_ 



1 



(3) 



taking into account that r+ = r_e^^^ ^°"^ofco(a+(i)^ rj.^^ 
complex reflection coefficient r_ = ^JRle^^ is computed 
numerically using the reflectivity R- of the left semi- 
lattice, given by the ratio of the reflected to the incident 
power, and the phase (\) of the wave reflected at the origin 
of the lattice. On the one hand, R- has to be used to 
prevent strong fluctuations in the local field £^^(0) due 
to the random distribution of the atoms arbitrarily close 
to the origin; on the other hand, the phase (\) needs to 
be taken at the origin of the lattice to avoid extra phase- 
shifts that appear at a distance in Gaussian beams (see 
the discussion below). As can be seen in Fig.gc), aPBG 
is observed in our model, which confirms the ability of 
this microscopic theory to capture the photonic structure 
of the atomic cloud. However, the LDOS is only in fair 
agreement with the results provided by the TM theory. 
The discrepancy in the phase ^ of the coefficient r_ of 
up to 7r/4 [see Fig. [4]^d)] is explained by the fact that 
our model treats the incident light as a Gaussian beam 
with a finite waist and space-dependent phase shifts (e.g., 
the Gouy phase), while the TM approach intrinsically 
assumes an incident plane wave. This affirms that Eq. ([3| 
is reliable only when the transverse finite-size effects are 
negligible. 




FIG. 4: (color online) Spectra of (a) the reflection coefficient 
i?, (b) the spontaneous emission, (c) the LDOS and (d) the 
phase of the complex reflection coefficient of the first half 
of the lattice. The simulations are realized for a semi-lattice 
with Nd — 50 and 150 layers of thickness a — 0.04Ao and 
the atomic density p — 5.95/Ao. The TM results correspond 
to the dashed black {Nd — 50) and the plain red (Nd — 150) 
lines, whereas the symbols stand for the simulations of the mi- 
croscopic model (black crosses and red circles, respectively). 
The latter have been done for the atomic layers of radius 9Ao 
filled with N — 3030 and 9090 atoms, respectively. 



Fig. |4] also shows that the band gap appears as the 
number of disks is increased, and larger systems may ex- 
hibit deeper band gaps in their spectrum. The simula- 
tion of large systems is actually the main limitation of 
our model since the complexity of the problem scales as 
. However, three-dimensional optical lattices typically 
contain several 10^ atoms [21], rendering the computa- 
tion feasible with standard computers. This makes our 
microscopic model particularly promising in the quest of 
three-dimensional omnidirectional photonic band gaps in 
optical lattices [151 Ull • 

Finally, we have verified that under some idealizing 
hypotheses the microscopic model boils down to the TM 
formalism which is commonly used to describe light prop- 
agation in one- dimensional atomic samples and charac- 
terize one-dimensional PBGs [9^. The first step involves 
a coarse-graining of the atomic structure describing the 
atoms as a continuous density distribution p(r). The 
atomic cloud is then characterized by a local refractive 
index n(r) = ^J^^^^~^K~pJ^^(lK^f^^^^ ^ and the wave 
propagating in it can be shown to satisfy the Helmholtz 
equation [22] : 



{y)\E = 



(4) 



Furthermore, assuming the system to be one- 
dimensional, which in practice means that its transverse 
size is much larger than the lattice period and the 
wavelength of the incident light, the scattering problem 
can be reduced to a one-dimensional wave-propagation 
problem and solved using classical techniques such as the 
TM formalism [23]. This explains the good agreement 
of the latter approach with our microscopic theory up 
to the point, where finite-size effects start playing a 
significant role [see Fig. [4]. According to our derivation, 
the reflection and transmission coefficients r^^ and 
for an atomic structure consisting of parallel disks of 
uniform density can be written in terms of the reflection 
and transmission coefficient amplitudes r and t for a 
single layer and the Bloch phase (pB'- 



r sin Nd(l)B 



sinTV^^B -tsm{Nd - 1)0b 
t sin (j)B 



(5) 
(6) 



sin Nd(l)B -tsin{Nd - 1)(/)b 
where 0b is determined by 

COS 05 = cos{kozd) cos{kza) — ^ — sin{kozd) sin{kza) 



with koz = ko cos and k^ = ko ^ v? — sin^ ^o- The 
detailed derivation of these results which are consistent 
with other models [9l [24j will be reported elsewhere. 

In conclusion, we have proposed a microscopic descrip- 
tion of the scattering of light from optical lattices and 
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demonstrated how multiple reflections from adjacent lat- 
tice sites can open a photonic band gap. The reconsider- 
ation of photonic bands under the microscopic scattering 
perspective leads not only to a deeper understanding of 
the phenomenon, but also offers a practical advantage 
of a larger range of applications than other models. For 
instance, our model naturally includes the cloud's granu- 
larity. Defects, such as site vacancies or flnite-size effects, 
can be easily taken into account, making the microscopic 
model particularly promising to study the transition from 
ordered to disordered clouds. This feature distinguishes 
our model from approaches based on the expansion of 
the electric fleld in terms of Bloch wave vectors and the 
solution of the Helmholtz equation in reciprocal space 
[151 US]- Those models most likely require optical lat- 
tices to be inflnite or perfectly periodic. 

Moreover, in comparison to the transfer matrix ap- 
proaches [9], the microscopic model can be readily ex- 
tended to two- or three-dimensional lattices of arbitrary 
geometry. It also allows for the description of exper- 
imental side effects in one-dimensional lattices, such as 
those considered in this paper, e.g., walk-off losses |lQ[ ll4j 
as the result of the flnite radial extension of the atomic 
disks, the deviation of the probe laser beam entering the 
atomic cloud caused by the refraction and the inhomoge- 
neous Stark-shift due to the intensity distribution of the 
trapping light [12]. 

The price to pay is that numerical simulations get 
quite heavy beyond a few 10^ atoms, which however, is 
not too far from experimentally relevant systems. These 
atomic numbers are proven to be sufficient enough to 
reach Bragg scattering and important reflection coefli- 
cients comparable to those obtained experimentally [13]. 
Hence, we believe the microscopic model to be appro- 
priate to characterize other collective phenomena usu- 
ally approached using coarse-grained theories, where the 
medium is described by a refractive index. In the case 
of disordered systems, the proposed model paves the way 
for a microscopic discussion of the extinction theorem [19] 
and the Abraham-Minkowski controversy [25 . 

Finally, it is worth noticing that, despite the fact that 
multiple scattering is naturally included into microscopic 
collective scattering models, the major part of recent 
studies avoid this regime, where the physical interpre- 
tation of the observed effects can be ambiguous. In con- 
trast, our work points out that ordered lattices repre- 
sent systems, where multiple scattering leads to the rel- 
atively simple and well-known phenomenon of photonic 
band gaps. 
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